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ABSTRACT 



In a previous study published in Astrobiology, we focused on the evolution of 
habitability of a 10 M® super-Earth planet orbiting a star akin to the Sun. This 
study was based on a concept of planetary habitability in accordance to the in- 
tegrated system approach that describes the photosynthetic biomass production 
taking into account a variety of climatological, biogeochemical, and geodynami- 
C^^ . cal processes. In the present study, we pursue a significant augmentation of our 

Q . previous work by considering stars with zero-age main sequence masses between 

0.5 and 2.0 Mq with special emphasis on models of 0.8, 0.9, 1.2 and 1.5 Mq. Our 
models of habitability consider again geodynamical processes during the main- 
^^ . sequence stage of these stars as well as during their red giant branch evolution. 

Pertaining to the different types of stars, we identify so-called photosynthesis- 
sustaining habitable zones (pHZ) determined by the limits of biological produc- 
tivity on the planetary surface. We obtain various sets of solutions consistent 
with the principal possibility of life. Considering that stars of relatively high 
masses depart from the main-sequence much earlier than low-mass stars, it is 
found that the biospheric life-span of super-Earth planets of stars with masses 
above approximately 1.5 Mq is always limited by the increase in stellar luminos- 
ity. However, for stars with masses below 0.9 Mq, the life-span of super-Earths is 
solely determined by the geodynamic time-scale. For central star masses between 
0.9 and 1.5 Mq, the possibility of life in the framework of our models depends 
on the relative continental area of the super-Earth planet. 
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Subject headings: extrasolar planets, geodynamics, habitable zone, planetary 
climate, stellar evolution and super-Earths. 



Introduction 



A central topic of contemporaneous astrobiology is the study of habitability around 
different types of stars, particularly main-sequence stars. Th is is motivated by the fac t that 



stars spent most of their lifetimes on the main-sequence (e.g.. lMaeder fc Meynetlll988l ). Fur- 
thermore, it is noteworthy that the distribution of stellar spectral types is strongly tilted 
toward low- mass stars, i.e., K and M-type stars, as implied by det ailed determinations of 
the in itial mass function in the solar neighborhood and beyond (e.g.. lKroupall2002l : IChabrier 
20031 ). For the present study we want to consider stars of large abundance in the galac- 
tic disk, within a reasonably large mass range, and which all have the potential of being 
central objects to habitable planetary systems. Thus, it is the aim of the present paper 
to significantly augment our previous study th at dealt with th e evolu tion of circumstellar 
habitability for stars akin to the Sun (1 Mq), see lvon Bloh et al\ (120091 ) . by considering stars 
with masses between 0.5 and 2.0 Mq. Moreover, we will adopt a wider grid of planetary 
masses, ecompassing 1 to 10 Earth masses (M®). 



von Bloh et all (120091 ) in Paper I presented a detailed thermal evolution model for a 10 



Earth-mass planet in the environment of a star like the Sun while also adopting an up-to-date 
solar evolutionary model (jSchroder fc Smithll2008l ). The planetary model of habitability has 
been based on the integrated system approach, which describes the photosynthetic biomass 
production taking into account a variety of climatological, biogeochemical, and geodynamical 
processes. This allowed us to identify a so-called photosynthesis-sustaining habitable zone 
(pHZ) determined by the limits of biological productivity on the planetary surface. The 
model also considered the principle possibility of habitability during stellar evolution along 
the Red Giant Branch (RGB). It was found that the solar pHZ increases in width over time 
and moves outward, as expected. For example, for ages of 11.0, 11.5, 12.0, and 12.1 Gyr, 
the pHZ is found to extend from 1.41 to 2.60, 1.58 to 2.60, 4.03 to 6.03, and 6.35 to 9.35 
AU, respectively. 

The approach of evaluating habitability based o n assessing the pHZ i s alternative to the 
study of the climatological habitable zone used by iKasting et all (119931 ) and others. Note 
that concerning Earth-mass planets, a detailed study of g eodynamic habitability based on 
the pHZ was previously presented by lFranck et all (l2000bl ). They found that Earth will be 
rendered uninhabitable after 6.5 Gyr as a result of plate tectonics, notably the growth of 
the continental area (enhanced loss of atmospheric CO2 by the increased weathering surface) 
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and the dwindling spreading rate (diminishing CO2 output from the sohd Earth). This work 
aheady considered systems of different types of main-se quence stars; Iiowever , it was based 
on an earher version of stellar evolution models given by lSchaller et all ( 119921 ) . Moreover, it 
also did not include stellar post-main-sequence evolution. 



The study by iFranck et al\ (J2000bl ) implies that there is no merit in investigating the 
future habitability of Earth during long-term stellar evolution, as in the framework of pHZ 
models, the lifetime of habitability is limited by terrestrial geodynamic processes. How- 
ever, this situation is expected to be significantly different for super-Earth pla nets due to 



inherent differences compared to Earth-mass planets (e.g.. [Valencia et a/.ll2007l ). A further 
motivation for this type of work stems from the ongoing discovery of su per-Earths in the 
solar neighborho od with the Gliese 876 ( iRivera et a/.ll2005l ) and Gliese 581 (jUdry et a/.ll2007l : 
Vogt et all I2OIOI ) systems as prime examples. In the following, we describe aspects of stel- 
lar evolutionary processes pertaining to our study. Next, we discuss the definition of the 
photosynthesis-sustaining habitable zone, including the relevant geodynamic assumptions. 
Thereafter, we present our results and discussion. Finally, we convey our summary and 
conclusions. 



2. Stellar Evolution Computations 



2.1. Methods 



As in Paper I, we base the selected evolution models on the well-tested Eggleton code, 
which allows us to follow the changes of stellar properties through the stages of the main- 
sequence, the red giant branch, and beyond. Our computations have been made with an 
advanced version of the Eggleton c ode, which c onsid e rs up dated opacities and an improved 
equation of state as described by iPols et al\ (1 19951 . Il998l ). For the abundance of heavy 
elements, which decisively affect the opacities, we use the near-solar value oi Z = 0.02. 
This choice is an appropriate representation of present-day samples of stars in the thin 
galactic disk, noting that these stars show a relatively narrow distribution {Z = 0.01 to 0.03) 
regarding heavy element abundances about this value. Besides other desirable characteristics, 
the adopted evolution code uses a self-adapting mesh and has a treatment of "overshooting" 
that has thoroughly been tested. Its two parameters, i.e., the mixing length and the overshoot 
length, have bee n cahbrated by utiliz ing giant and supergiant stars in well-studied, eclipsing 
binary systems (ISchroder et a/.lll997l ). 



The evolution code, as used by us, also considers a detailed description of the stellar 
mass loss, including its impact on the evolution during the stages of giant star evolution 
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foUowing ISchroder fc Cuntzl (J2005[ ). The treatment of the stellar mass loss rate as function 
of the governing stellar param eters has successfully bee n tested on globular clusters as well as 
for a set of well studied stars (ISchroder fc Cunta 120071 ). Thus, the attained models of stellar 
evolution are expected to provide an accurate description of the time-dependent behavior of 
stellar luminosity (see Fig. 1), as well as other stellar quantities, from the main-sequence to 
the RGB stage, for the entire range of masses considered in our study. Also note that as a 
consequence of the steadily increasing mass loss, particularly on the upper RGB of various 
stars including the Sun, the orbital distance R of any putative planet increases as i? oc M~^ 
with M^K as stellar mass, owing to the conservation of the orbital angular momentum of the 
planet. 



2.2. Mass range of interest 



In Table 1, we present some characteristics of stellar evolution models for the mass 
range of 0.5 to 2 Mq. F or comparison, we note th at the solar age, according to our models, 
is given as 4.58 Gyr, see ISchroder fc SmithI (120081 ) (for a metallicity oi Z = 0.018), and the 
Sun's present-day effective temperature is identified as T^s = 5774 K. Masses below 0.5 Mq, 
i.e., most M-type stars, have been omitted from our study. Their stellar evolution can be 
neglected, luminosity and stellar temperature remain unchanged. 

Moreover, stars of relatively large masses, i.e., M > 2 Mq are also considered poor 
candidates for supporting life because the stellar lifetimes are almost certainly too short to 
permit the onset of biology. In Table 1, we list information on the target stars taken into 
consideration. Particularly, we convey the time for each stellar model at which the initial 
(i.e., zero-age main-sequence) luminosity I/zams has increased by 50%, and when Lzams has 
doubled. Stars above 2 Mq reach the stage at which they leave the main-sequence phase 
after less than a billion years. If we set a relatively low age limit of, say, 300 million years as a 
requirement for the development of complex life forms during the phase of slowly increasing 
stellar luminosity, we would limit the mass range of potential host stars of life-bearing planets 
to a maximum of about 1.3 Mq. 

In Table 1, we also list the stellar effective temperatures, which are a steep monotonic 
function of the stellar mass. Note that due to the appearance of increased photospheric UV 
fluxes, relatively high effective temperatures by th emselves repr e sent an adverse fa ctor for 



the general possibility of advanced life forms (e.g., lCockelllll999l : ICuntz et a/.ll2010l ). which 
is particularly relevant for stars of spectral type mid-F and earlier, although the biological 
impact of stellar UV environments is typically significantly altered by the attenuation of 
planetary atmospheres. This consideration is a further motivation to disregard stars with 
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masses beyond 1.5 or 2.0 M© (i.e., early F stars) in the context of exobiology. Towards the 
low mass limit at about 0.5 M©, there are virtually no stellar evolutionary changes for these 
stars while being on the main-sequence owing to the current age of the Universe. However, 
regarding M-type stars, adverse influences on the origin and development of sophisticated 
life forms also exis t , particularly associated with strong UV, EUV and X -ray flaring (e.g., 
Hawlev et a/lboosl : IScalo et alhoQli ICuntz et a/.lboiol : ISeeura et a/.ll201oh . 



Close to about one solar mass, we find an interesting result: Stars of, e.g., 0.9 Mq, still 
spend many billion years in increasing their luminosity from 1.5 times to 2 times of their zero- 
age main-sequence luminosity. In that phase, which is sufficiently stable and long-lived to 
provide the development of life forms, the initially significantly cooler effective temperature 
of such a star is closing in on that of the early Sun. This change of quality and quantity of the 
host star radiation could "unfreeze" a planet, which was initially outside the habitable zone, 
and still leave it with enough time to develop complex life forms. When reaching about 1.5 
times their zero-age main-sequence luminosity, stars of 0.9 solar masses have reached an age 
of about 8 billion years (see Fig. 2), much like the older stars of the thin galactic disk. Hence, 
several such cases should be found in present-day stellar samples of the solar neighbourhood. 
We also note that these same stars continue to evolve quite slowly. It will take them several 
further billions of years to leave the main-sequence and to advance on the RGB. 



2.3. A closer look at the mass range of 0.8 to 1.5 Mq 

As previously pointed out, the most promising types of stars for hosting life-bearing 
planets regarding our grid of models are those of 0.8, 0.9, 1.0 (Sun), 1.2 and 1.5 Mq. All these 
stars have a lot in common: they all evolve into red giants (RGB phase) after central hydrogen 
burning has ceased. The expansion of the outer layers of a RGB star is driven by two related 
processes: the contraction of the inner helium core, which has lost its stabilising energy 
production, and the growing energy production by the hydrogen burning shell around it. This 
layer right above the contracting core experiences a steady density increase, which drives up 
the resulting stellar luminosity. The now much higher energy output of the hydrogen burning 
shell as well as the much higher density of the core region lead to a very expanded outer 
stellar structure, associated with a significant increase in stellar luminosity and radius (see 
Fig. 2). Additionally, the stella r structure is shaped by the onset of a cool wind, associated 
with significant mass loss (e.g.. iDupree fc Reimerd 119871 ). 



For our models of 1.0, 1.2 and 1.5 Mq, the appearance of a helium fiash marks the end 
of the RGB phase, as the stars settle into a new equilibrium with a much less compact core 
(due to its new source of energy) and a, consequently, much less expanded outer structure. 
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However, the post-RGB stellar evolution does not offer any promising aspects for the origin 
of life, since the relatively stable phase of central helium burning lasts only a few hundred 
million years. Furthermore, any ancient life forms on previously habitable worlds will most 
likely be destroyed, since in the extreme RGB phases of those stars amount to 200-300 times 
its initial size and several thousand times its initial luminosity. This entails a complete 
nullification of previous zones of circumstellar habitability established during stellar main- 
sequence evolution. 

However, despite such remarkable homogeneity concerning their main physical processes, 
stars in the range of 0.8 to 1.5 solar masses differ largely in their lifetimes; see Fig. 2 and 
Table 1 for details. On the upper RGB, the much longer times spent in that phase by 
low-mass stars have an interesting consequence: their mass loss experienced in this phase 
amounts to more than 0.4 Mq. This has dramatic implications for any existing planets at his 
stages noting that their orbital distance R increases as i? oc M~^ . This may imply that an at 
that time habitable planet may permanently leave the stellar habitable zone or, conversely, 
an i nhabitable p l anet m ay enter the stellar habitable zone, and becomes "unfrozen"; see. 



Lopez et al\ (120051 ) and Ivon Bloh et al\ (120091 ) for previous results for stars akin to the 



e.g., 
Sun. 

In summary, it is the stellar main-sequence phase that is most relevant for the evolution 
of life on planets at a suitable distance. The slower the star increases its luminosity, the 
longer a planet stays within the slowly expanding habitable zone. Thus, all stars between 
0.9 and 1.2 Mq masses are potential host stars for life from the beginning of stellar evolution. 
Stars between 0.8 and 0.9 solar masses show an even slower change in stellar effective tem- 
perature and luminosity; in fact, they reach solar luminosity after 8 billion years. Detailed 
studies for the evolution of habitability, based on the def inition of the clima t ologic al habitable 
zones (see below), durin g stellar lifetimes were given by lUnderwood et al\ (120031 ) and others. 
Jones et al\ (120051 . |2006| ) applied this framework to known exoplanetary systems to evaluate 
the possibility of habitable "Earths" . Stella r evolutionary aspects are a lso relevant to the 
future of life on Earth, as discussed by, e.g., ISackmann et al\ (Il993l ) and lSchroder fc Smith 
( 120081 ) ■ although the impact of geody namic processes app ears to be dominant for setting the 
time limit of terrestrial habitability ( JFranck et a/.ll2000bl ). 
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3. Habitability of Super-Earth Planets 

3.1. Photosynthesis-sustaining habitable zone (pHZ) 

Customary simulations about the role of stellar evoluti on on circumstellar habitability 
typic ally invoke the framework of the climatic habitable zone ( Kasting et a/JlQQSt ISelsis et al. 
20071 ). In this type of models, the zone of habitability at a given evolution ary time for a star 



with l uminosity L and effective temperature Teg can be calculated following lUnderwood et al. 
J2003h as 

L \5 „ f L 



Rh 



-^O ■ '-'ini-'c 



cffJ 



R, 



■out 



Lq ■ «Sout(-ic 



off J 



(1) 



with 5'in(Tcfj) and ^outl^efr) described as second order polynomials. 



It is the purpose of this paper to expand the focus of Paper I, which is solely focused 
on the Sun. Paper I is based on the concept that the habitability of super-Earth planets is 
evaluated via a modified Earth-system model that describes the evolution of the temperature 
and atmospheric CO2 concentration. On Earth, the carbonate-silicate cycle is the crucial 
element for a long-term homeostasis under increasing solar luminosity. On geological time- 
scales, the deeper parts of the Earth are considerable sinks and sources of carbon. 

Our numerical m odel was previously applied to the su per-Earths Gl 581c a nd Gl 58 Id 



(Ivon Bloh et a/.l 120071 ). the putative super-Earth Gl 58 Ig 



fictitious Earth-mass planets for 47 UMa and 55 Cnc (jCuntz et al 



von Bloh et al. 



20031 



20111) as well as 



Franck et al. 2003 



von Bloh et a/.ll2003l ). among numerous other studies. This model couples the stellar lu- 
minosity L, the silicate-rock weathering rate F^r and the global energy balance to obtain 
estimates of the partial pressure of atmospheric carbon dioxide -Pco2! the mean global surface 
temperature Tsurf, and the biological productivity 11 as a function of time t (Fig. [3]). The 
main point is the persistent balance between the CO2 sink in the atmosphere-ocean system 
and the metamorphic (plate tectonic) sources. This is expressed through the dimensionless 
quantities 

/wr(t)-/A(t) = /sr(t), (2) 

where fwrit) = -Fwr(^)/-Fwr,o is the weathering rate, fA(t) = Ac(t)/Acfi is the continental area, 
and fsrit) = S{t)/So is the areal spreading rate, which are all normalized by their present 
values of Ea rth. Eq. (|2]) can be rearranged by introducing the geophysical forcing ratio GFR 
(1Volklll987h as 

UriTsuruPco,) = ^ =: GFR(t) . (3) 

JA 

Here we assume that the weathering rate only depends on the global surface temperature 
and the atmospheric CO2 concentration. For the investigation of a super-Earth under ex- 



ternal forcing, we adopt a model planet with a prescribed continental area. The fraction of 
continental area with respect to the total planetary surface Ja is varied between 0.1 and 0.9. 

The connection between the stellar para meters and the planetary climate can be ob- 



tained by using a radiation balance equation (IWilliamslll998l ) 



AttR^ 



[1 - a(Ts,rf, PcoJ] = 4/R(Ts,rf, PcoJ, (4) 



where a denotes the planetary albedo, //? the outgoing infrared flux, and R the distance from 
the central star. The Eqs. ([3]) and (jl]) constitute a set of two coupled equations with two 
unknowns, Tguri and Pco2! if the parameterization of the weathering rate, the luminosity, 
the distance to the central star and the geophysical forcing ratio are specified. Therefore, a 
numerical solution can be attained in a straightforward manner. 

Resulting from the described model the domain of distances R for a super-Earth planet 
from its respective main-sequence star for which a photosynthesis-active biosphere is pro- 
ductive (n > 0) can be calculated. It is defined as the photosynthetis-sustaining habitable 
zone given as 

pHZ(t) := {R\U{Pco,{R,t),T,^,i{R,t)) > 0}. (5) 

In our model, biological productivity is considered to be solely a function of the surface 
temperature and the CO2 partial pressure in the atmosphere. Our pa rameterization yield s 



zero productivity for T^^,f < 0°C or T^,,,f > 100°C or Pcoa < lO^^ bar flFranck et a/.ll2000af ). 
The inner and outer boundaries of the pHZ do not depend on the detailed parameterization 
of the biological productivity within the temperature and pressure tolerance window. 

The maximum life span of a photosynthetic-active biosphere tmax is the point in time 
when a planet at an optimum distance from its central star i?opt finally leaves the pHZ 

Wx = max|pHZ(t)| >0. (6) 

Note that the evaluation of pHZ(t) and tmax will be essential to our models of circumstellar 
habitability. 



3.2. Comments on the thermal evolution model 

Parameterized convection models are the simplest models for the investigation of the 
thermal evolution of terrestrial planets and satellites. They have su ccessfully been applied 



to the evol ut ion of Mercury, Venus, E arth, Mars, and the Moon (IStevenson et all Il983 



SleepI |2000[ ) . iFranck &: Bounamal (Il995[ ) studied the thermal and volatile history of Earth 
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and Venus in the framework of comparative planetology. The internal struc ture of massive 



terres trial planets with one to ten Earth masses has been investigated by [Valencia et al 



( I2OO6I ) to obtain scaling laws for total radius, mantle thickness, core size, and average density 
as a function of mass. Similar scaling laws were found for different compositions. We will 
use these scaling laws for mass-dependent properti es of super-Earths and als o the mass- 
independent material properties previously given by iFranck fc Bounamal ( 1l995l ). 



The thermal history and future of a super-Earth planet has to be determined to cal- 
culate the spreading rate for solving key Eq. ([2]). A parameterized model of whole man- 
tle convection including the volatile ex change between the mantle and surface reservoirs 
( IFranck fc Bounamal Il995l : iFranckl Il998[ ) is applied. The key equations used in our present 
stu dy are in accor d with our p revious work focused on Gl 581c, Gl 581d, and Gl 581g; 



see 



von Bloh et al\ ( 120071 . 1201 if ) for details. A key element is the computation of the areal 



spreading rate S] note that S" is a function of the average mantle temperature T^, the sur- 
face temperature Tsnrfi the heat flow from the mantle qm-, and the area of ocean basins Aq 
( jTurcotte fc Schubertll2002l ). It is given as 



S 



q^TTKAo 



Ak\T^-T, 



surf J 



\2 ' 



(7) 



where k is the thermal diffusivity and k the thermal conductivity. To calculate the spreading 
rate, the thermal evolution of the mantle has be to computed: 



^vrpc(i?^-i?3)^ 



-A7,Rlqm + -'KE{t){Rl- R, 



% 



where p is the density, c is the specific heat at constant pressure, E is the energy production 
rate by decay of radiogenic heat sources in the mantle per unit volume, and Rm and Re are 
the outer and inner radii of the mantle, respectively. To calculate the thermal evolution for a 
planet with several Earth masses the planetary parameters have to be adjusted. Therefore, 
we assume 



Rp 

Rff 



M 



0.27 



(9) 



and where Rp is the planetary radius, see [Valencia et all (120061 ). 



The total radius, mantle thickness, core size and average density are all functions of 
mass, with the subscript © denoting Earth values. The exponent of 0.27 has been obtained 
for super-Earths. The values of Rp, Rm, Re, as well as the other planetary properties are 
scaled accordingly. Table 2 gives a summary of t hese size pararn e ters f o r the planetary 
models between 1 and 10 M®; see the studies by Ivon Bloh et all (120071 . I2OIII ) for addi- 
tional information including parameters for the mass-independent quantities. According to 
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Valencia &: O'ConnelJ ( 120091 ). we assume that a more massive planet is likely to convect in 
a plate tectonic regime similar to Earth. Thus, the more massive the planet, the higher the 
Rayleigh number that controls convection, the thinner the top boundary layer (lithosphere) , 
and the higher the convective velocities. 

We recognize that there is an ongoing deba te about the onset and significance of plate 
tecton ics on super-Earth planets. On one hand, IStein et al\ (120041 ) and lO'Neill fc Lenardic 
( 120071 ) pointed out that there might be in an episodic or stagnant lid regime, considering 
that increasing the planetary radius acts to decrease the ratio of driving to resisting stresses, 
an effect that appears to be robust when increases in planetary grav ity are included. Other 



thermodynamic processes that peak at different pl a net rn asses (e.g.. iNoack fc Breuer 



20111 



1) 



operate as well. On the other hand, [Valencia et al\ (120071 ) and [Valencia fc O'Connelll (J2009) 
argued that planet tectonics of super-Earths should be considered inevitable. Their models 
show that as the planetary mass increases, the shear stress to overcome resistence to plate 
thickness increases while the plate thickness decreases, thereby enhancing plate weakness. 
These effects contribute favorably to the subduction of the lithosphere (like on Earth), a 
crucial component of plate tectonics. Noting that in the framework of our model, a plate- 
tectonic-driven carbon cycle is considered necessary and essential for carbon-based life, we 
will assume the existence of plate tectonics on super-Earths for our models in the following. 



4. Results and Discussion 



4.1. Habitability based on the integrated system approach 

We calculated the photosynthesis-sustaining habitable zones for a 10 M^ super-Earth 
with relative continental areas between r = 0.1 and 0.9 orbiting central stars in the mass 
range between 0.8 and 1.5 Mq. The results are shown in Fig. HJ The width of the pHZ 
during the main-sequence evolution is found to be approximately constant for all stellar 
models taken into account. For stellar mass M < 0.9 Mq the pHZ ceases to exist before the 
main-sequence evolution ends, while for higher stellar masses the width of the pHZ increases 
significantly and moves outwardly in time. The time when this processes occurs critically 
depends on the stellar mass, as expected. 

"Water worlds", i.e., planets mostly covered by oceans are favored in the facilitation 
of habitability during this stage as previously obtained in the context of long-term evo- 
l ution for hypothet i cal Earth-mass and super-Earth planets in various star-planet systems 
(JFranck et a/.ll2003l : Ivon Bloh et a/.ll2003l . 120071 . l2009l ). In principle, the favored existence of 
water worlds instead of land worlds around stars during their phase of red giant branch evo- 
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lution constitutes a possible extension of the outer edge of circumstellar habitability. The 
reason for the favorism of water worlds under those circumstances is that planets with a 
considerable continental area have relatively high weathering rates that provide the main 
sink of atmospheric CO2. Therefore, this type of planets are unable to build up C02-rich 
atmospheres that prevent the planet from freezing, which in turn is thwarting the possibility 
of photosynthesis-based life. 

The maximum life spans tmax (Eq- ED of super-Earth planets have been calculated for 
a grid of planetary masses between 1 and 10 M^ with increments of 1 M^ and relative 
continental areas r between 0.1 and 0.9. Figure ([5]) depicts tmax as a function of the planetary 
mass. It is found that tmax increases with planetary mass and, furthermore, decreases for 
increasing r values. In addition, we depict the time-spans as a function of the stellar mass 
when the central star reaches a luminosity of L = 1.5 Liams and L = 2 Liams with Lzams 
as the initial (i.e., zero-age main-sequence) luminosity at time t = 0. 

Moreover, in Fig. [6] the critical central stellar mass is depicted dependent on the plan- 
etary mass up to which the life span of the biosphere is solely determined by tmax (colored 
shaded areas) instead of the increase in stellar luminosity (white area) . It is found that the 
biospheric life span of super-Earth planets for central stars with masses above about l.SM© 
is always limited by the increase in stellar luminosity. For central star masses below O.QM© 
it is solely determined by tmax, i-e., the maximal geodynamic life span of the biosphere can 
always be realised irrespectively of the stellar nuclear evolution. For central star masses 
between 0.9 and 1.5Mq the situation depends on the relative continental area r. 

Main-sequence stars with masses between 0.5 and 0.9 Mq belong to late-type G, K, 
and early M stars. The luminosity of such stars remains almost constant during their en- 
tire nuclear evolution on the main- sequence. Therefore, planets within the habitable zone 
would remain habitable for a very long period of time if the geodynamic evolution of the 
planet is not taken into account. The planetary evolution is driven by the long-term cooling 
of the planetary interior. Therefore, diminished internal forcing will affect the planetary 
habitability of the biospheres. In the framework of the adopted models it is found that the 
maximum life span of super-Earth planets around low mass stars depends almost entirely 
on the properties of the planet. It is known, however, that other factors are also expected 
to potentially impact or limit the habitability in the environments of late K and M stars, 
encompassing effects associated with, e.g., tidal locking of the planet and stellar flares (e.g. 



Lammei 



Heller et al. 



2OIII). 



20071: iTarter et a/.l 120071 : ISelsis et a/.l 120071 : iLammer et q/.II2009I : ISegura et a/.ll2010 
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5. Summary and Conclusions 



We investigated the habitability of super-Eartli planets in the environments of main- 
sequence stars of masses between 0.5 and and 2.0 M©. This work c omplements a previou s 



study that has been focused on star akin to the Sun (M = IM©), see lvon Bloh et al\ ( 120091 ). 
Both studies employ the concept of planetary habitability based on the integrated system 
approach that describes the photosynthetic biomass production taking into account a variety 
of climatological, biogeochemical, and geodynamical processes. Among the various assump- 
tions entertained as part of our study, we consider that super-Earth planets are expected to 
develop plate tectonics, an assumption critical to our models, although it is currently under 
scrutiny. Pertaining to the different types of stars, we identify so-called photosynthesis- 
sustaining habitable zones (pHZ) determined by the limits of biological productivity on the 
planetary surface. We obtain various sets of solutions consistent with the principal possibil- 
ity of life. Our models are based on advanced version of the Eggleton stellar evol ution code 



that conside rs updated opacities and an updated equation of state as described by lPols et al 



( 119951 . Il998l ) . Additionally, an i mpro ved description of mass loss is implemented following 



the work by ISchroder fc Cuntzl (120051 ) . 



Considering that high-mass stars depart from the main-sequence much faster than low- 
mass stars, it was found that the biospheric life-span of super-Earth planets of stars with 
masses above approximately 1.5 Mq is always limited by the increase in stellar luminosity 
dictated by stellar nuclear evolution. However, for stars with masses below 0.9 Mq, the 
life-span of super-Earths it is solely determined by the geodynamic time-scale. Clearly, the 
luminosity of stars with spectral type mid-K or later has not deviated much from their ZAMS 
luminosities owing to the limited age of the Universe. For central star masses between 0.9 
and 1.5 M©, the possibility of life in the framework of our models depends on the relative 
continental area of the super-Earth planets. 

A crucial finding of our study is that planets of large geological age are expected to 
be water worlds, i.e., planets with a relatively small continental area. The reason why 
land worlds are considered to be unfavorable is that planets with a considerable continental 
area have relatively high weathering rates that provide the main sink of atmospheric CO2. 
Therefore, this type of planets are unable to build up C02-rich atmospheres that prevent the 
planet from freezing or allowing for the possibility of photosynthesis-based lif e. Our study 



i s also both expanding and superseding previous work, including the study by iFranck et al 



( l2000bl ). This work solely considered Earth-mass planets instead of also including super- 
Earth planets. Additionally, it considered a linear growth model for the continental growth 
and has refrained from taking into account post-main sequence evolution. Although the 
evolutionary time-scales for stars on the RGB are short compared to the time spent on 
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the main- sequence, habitable planets ar ound red gia n ts sho u ld be considere d a re alistic 



possibility, as previously pointed out by iLopez et al\ (120051 ) . I von Bloh et al\ (120091 ) . and 



others. However, for central stars with higher masses than the Sun, a more rapid evolution 
will occur that will also place the significant temporal and spatial constraints on planetary 
habitability when the central stars have reached the RGB. 
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Fig. 1. — Hertzsprung- Russell diagram (logL versus logTgff) of stellar evolution tracks with 
initial masses of 0.8, 0.9, 1.0 (Sun), 1.2 and 1.5 M©, reaching from the ZAMS to the tip 
of the RGB. In the first two cases, the stars lose so much mass on the RGB that they do 
not ignite He-burning but directly evolve into white dwarfs. The AGB branches of the more 
massive stars have been omitted for the sake of clarity. 
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Fig. 2. — Stellar evolution models following ISchroder fc Smithl ( l2008l ). depicting the lumi- 
nosity (solid line), the effective temperature (dashed line), and the mass (dash-dotted line). 
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Fig. 3. — Box model of the integrated system approach adopted in our model. The arrows 
indicate the different types of forcing, which are the main feedback loop for stabilizing the 
climate (thick solid arrows), the feedback loop within the system (thin solid arrows), and 
the external and internal forcings (dashed arrows). 
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Fig. 4.— The pHZ for a 10 Earth-mass planet orbiting an M = 0.8, 0.9, 1.0, 1.2, 1.5 M© 
star also considering post-main sequence evolution. The gray shaded area at the bottom 
denotes the time period of the stellar main-sequence evolution. Note the "bending-over" of 
the color-coded areas for stars with masses of 1.0 Mq and above due to the drastic increases 
in the stellar luminosities. 
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Fig. 5. — The maximum life span of the biosphere tmax in dependence on the planetary mass 
M (colored solid lines). The graphs are color coded by the assumed portion of the planetary 
surface covered by continents r. The solid black line denotes the time when the star reaches 
a luminosity of L = 1.5 Lzams, while the dashed hne denotes the time for L = 2.0 Lzams- 
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Fig. 6. — The critical mass of the central star as a function of the planetary mass up to 
which the life span of the biosphere is solely determined by the maximum geodynamic life 
span (color-coded area) instead of the time of the stellar nuclear evolution, i.e., when the star 
reaches a luminosity of L = 2.0 Lzams (white area). Different colors correlate to different 
relative continental areas of the model planet. 
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Table 1. Stellar Parameters 



M 


log(LzAMs)/^0) 


TeS 


Age(1.5 X LzAMs) 


^cff 


Age (2 X LzAMs) 


^cff 


(Mo) 




(K) 


(Gyr) 


(K) 


(Gyr) 


(K) 


2.0 


1.204 


9160 


0.98 


7390 


1.20 


7100 


1.9 


1.113 


8790 


1.18 


6975 


1.39 


6885 


1.8 


1.017 


8400 


1.45 


6455 


1.63 


6380 


1.7 


0.914 


7990 


1.84 


6130 


1.92 


6260 


1.6 


0.802 


7555 


2.15 


6080 


2.28 


6340 


1.5 


0.679 


7100 


2.6 


6000 


2.8 


6300 


1.4 


0.543 


6710 


2.6 


6220 


3.0 


6180 


1.3 


0.391 


6420 


3.3 


6145 


3.7 


6165 


1.2 


0.221 


6140 


3.6 


6110 


4.7 


6070 


1.1 


0.015 


5820 


3.8 


5990 


6.1 


5960 


1.0 


-0.163 


5580 


5.5 


5750 


8.8 


5770 


0.9 


-0.383 


5180 


8.0 


5430 


12.3 


5530 


0.8 


-0.625 


4735 


12.5 


5015 


18.5 


5185 


0.7 


-0.897 


4260 


21.0 


4540 


30.0 


4720 


0.6 


-1.170 


3930 


37.8 


4110 


53.1 


4260 


0.5 


-1.430 


3750 


72.5 


3865 


99.8 


3960 
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Table 2. Stellar and Planetary Parameters 



Parameter 




Value 




Unit 


Description 




IMe 


2 Me 


5 Me 


10 Me 






9 


1.00 


1.38 


2.10 


2.88 


^e 


gravitational acceleration 


R, 


6378 


7691 


9849 


11,876 


m 


planetary radius 


IhQ 


3471 


4185 


5360 


6463 


m 


inner radius of the mantle 


Rm 


6271 


7562 


9684 


11,677 


m 


outer radius of the mantle 



